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Abstract 

The production near threshold of isoscalar pion pairs in the pd — > 
3 He (tt 7r)° reaction is estimated in a two-step model which successfully 
describes the production of rj, oo and rj' mesons. A virtual pion beam, 
generated through an N N — > dir reaction on one of the nucleons in 
the deuteron, produces a second pion via a tt N — > tt tt N reaction on 
the other nucleon. Using the same scale factor as for heavy meson 
production, the model reproduces the total 7r°7T production rate de- 
termined at an excess energy of 37 MeV. There are some indications 
in the data for a suppression of events with low tttt masses, as in the 
7r~ p — > 7T° 7T° n reaction, and this is confirmed within the model. The 
model suggests that a significant fraction of the charged pion pro- 
duction in the pd — > 3 He 7r + 7r~ reaction at Q = 70 MeV might be 
associated with isoscalar pion pairs, though this does not explain the 
strong dependence observed on the tt + tt~ relative momentum angle. 
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The study of neutral two-pion production through the pd — > 3 HeX° 
reaction has a long history. At excess energies Q (the cm. kinetic energy 
in the final state) around 200-300 MeV, sharp structure is seen at missing 
masses of about 310 MeV/c 2 [1], [|. The absence of any significant strength 
in the pd — > 3 HX + channel at low mx means that the effect is associated 
with isospin-zero pion-pion pairs which, because of the available energy, must 
be dominantly in s-waves. Although a quantitative explanation of the ABC 
enhancement has not yet been provided for this reaction, a similar effect 
in np ^ dX° has been shown to originate from the excitation of two A- 
isobars || . The prominent ABC peaks in the dd — > 4 He X° case have also 
been shown to be due to double pion p-wave production ||. 

The experimental picture changes dramatically at lower energies. For 
values of Q around 70-90 MeV, the angular distributions observed by the 
MOMO group || U for the exclusive pd —>■ 3 He7r + 7r~ reaction suggest 
strongly that the ir + ir~ spectrum is mainly p-wave in nature, and hence has 
isospin-one. The measurement of the 7r°7r°/7r + 7r~ charge ratio at CELSIUS 
at an excess energy (with respect to the 7i°ir° threshold) of Q = 37 MeV |7j 
shows that there is significant 1 = 1 production even at this much smaller Q. 
The isoscalar 7r° tt° spectrum, determined either by direct measurement |7| or 
through the subtraction of exclusive 3 He7r + tt~ data from an inclusive pd — > 
3 HeV° measurement 0, shows that there is no s-wave ABC enhancement 
at low Q. On the contrary, there are rather indications that the s-wave cross 
section is actually suppressed at low 7r° tt° masses as compared to phase 
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space 0. It is the aim of the present paper to demonstrate that near- 
threshold isoscalar two-pion production in the pd — > 3 He ir n can be described 
in terms of sequential single-pion production. 




Figure 1: Dynamical model for the pd —>■ 3 He(7T7r)° reaction in terms of 
sequential pp — ► dir + and 7r + n — > (7T7r)°p processes. There is an analogous 
contribution from intermediate neutral pions. 

The large momentum transfers required to produce heavy mesons, such as 
the T] or u>, through the pd — > 3 HeX° reaction mean that two-step processes 
which minimise the momentum mismatch in the nuclear wave functions can 
provide the dominant driving force. In one such model, a pion produced 
on one of the nucleons in the target deuteron is converted into the observed 
heavy meson through an interaction on the second of the target nucleons |§ . 
Apart from an ad hoc overall normalisation factor N rs 2.4, which may reflect 
the retention of only bound intermediate deuteron states in the calculation, 
this approach describes well the threshold amplitudes for producing i], to 
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and T]', though the experimental <j) yield is a little too high ||. This success 
may be attributed to the fact that the intermediate pion in the diagram 
is close to its mass shell. We wish to apply the same model to isoscalar 
mi production by introducing rather a final tt N — > (tttt) N process, as in 
Fig. 1. For definiteness, we consider 7r°7r° production; estimates of charged 
pion production in the 1 = channel will then follow from isospin invariance, 
after correcting for the pion mass difference. 

The principal difference with the earlier work || is that the low mass 
7T° 7T° system is in a + state and so there is a parity change at the tt + n — > 
7T° 7T° p vertex. Parameterising this amplitude in terms of two-component 
Pauli spinors u p (n) as 

M(vr+n -> n°7r°p) = a(m n7T , W nN ) uj a • u n , (1) 

the corresponding differential cross section for s-wave production is 

da(n + n -> 7r°7r°p) = da(rr'p -> 7r°7r°r2) = — — - — T - |a(m^, W vN )\ 2 fc* dm nT , 

(2) 

Here p and p' are the incident and final nucleon momenta in the overall 
cm. system where the total energy is W n ^. In the 7T7T rest frame, k* is the 
relative momentum, which is related to the nir invariant mass through m wn = 
2^Jk* 2 + ml- We shall neglect the angular dependence of the amplitudes in 
the present work. 

Because of the nature of the two-step process in Fig. 1, only small Fermi 
momenta are required. Working to first order in these momenta, as in ||, 
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we find that the amplitudes are proportional to the complex form factors 

POO 

S a/3 {W,V) = (2tt) 3 / dte itASo C(-tW)^(tV) . (3) 

These involve integrals over configuration-space deuteron (cpp) and 3 He (ip a ) 
wave functions, where a, /3 — (0,2) represent nuclear S- and D-state com- 
ponents. The energy mismatch AE between the intermediate and external 
energies for zero Fermi momenta is generally small for near-threshold heavy 
meson production in this model. 

The relativistic relative velocity vectors V and W, 

1 



v = 2 1 Pmr 1 

3 ^(Iptttt — |pd) 2 
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w = 4 
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where p^ is the total tttt momentum vector and p^ that of the initial 
deuteron in the overall cm. frame. The component of V along pa must 
be subjected to a Lorentz contraction ||. The relativistic energies E^ are 
evaluated at the values of the momenta indicated. 

For zero Fermi momenta, the pp — > dir + amplitudes should be evaluated 
in the forward direction for threshold heavy meson production. The forward 
direction assumption is also very good even away from threshold provided 
that pa ^> p nn , as it is in cases under investigation. There are two pp — > dn + 
amplitudes in the forward direction but, at the energies required here, the 
helicity-zero completely dominates over the helicity-one Keeping then 



6 



only the dominant amplitude A, the cm. differential cross section is 

The evaluation of the unpolarised differential cross section in this model 
is similar to that for the production of single heavy mesons || and leads to 

x \A\ 2 \a(m^,W nN )\ 2 klpldm^dfl Ke , (6) 



where an isospin factor of | has been included to account for the 7r°-exchange 
term in Fig. 1. The form factor combinations required are 

S a = Sqo — S20 V2 , Sb = S02 — S22V2 . (7) 

In order to describe (rj, u, rf, <p) production, it was found necessary to mul- 
tiply the analogous prediction by a normalisation factor N = 2A |J. 

Single pion production in pion-nucleon collisions has been measured in 
many charge states near threshold and parameterisations given for the total 



cross sections as functions of the beam energy pT| . The charge dependence 
indicates that the cross section is dominated by I = pion pairs for Q < 
100 MeV. This is consistent with the smallness of the anisotropy in the 
angular distribution of the tt + tc~ relative momentum for ir~ p — > n + tt~ n, 
which arises from s-p, and hence 1 = 0/1=1 interference |12| . Data on 



the 7r p — ► 7r 7r n reaction in the Q ~ 50 — 100 MeV region show clear 



evidence for the suppression of events with low [13|. This is also seen 



for 7r p — » 7r + 7r n but not n + p ^ n + n + n, where only 1 = 2 pion pairs 



are produced jnj. In the dynamical model of the Valencia group this 
shift towards higher due to an accidental cancellation between two 

contributions, one of which involves the double pion p-wave decay of the 
Roper resonance iV*(1440) — > iV7r 7r . 



Although the production of isospin-two ixix pairs is small [ I2| , it has to 
be subtracted from the 7r 7r data to get the 1 = rate required in Eq. flU) 
This subtraction is model-dependent and, for this purpose, we have used 



the predictions of the Valencia model |TJ]], which describes reasonably well 



the shape of the experimental data (13], [T^]. A global fit to their 1 = 
predictions, renormalised slightly to agree with the overall 7r°7r° amplitude 



analysis of Lowe and Burkhardt flll ] , gives 



\a{m wn , Q)\ 2 = (1.092 - 0.0211Q + 0.00015Q 5 



64vr 3 

+ (4.18 + 0.0075Q - 0.00098Q 2 ):r + (47.65 - 0.935Q + 0.00743Q 2 ):r 2 /xb/MeV 2 , 

which is valid up to Q' = Q (1 +m 7r /m Hc ) 2 /(l + m n /m p ) 2 « 100 MeV. Given 
that x = m nn /m 7T — 2, this illustrates the suppression of the matrix element 
at low m^ w , a feature which becomes even more pronounced at higher Q. 

Our predictions for the pd — > 3 He(7T7r)° total cross section divided by 
Q 2 , obtained using the same value of iV = 2.4 which gave good agreement 
for heavy meson production, are to be found in Fig. 2. The steady increase 
with Q is mainly a reflection of the energy dependence of the pp —>■ d-K + and 
7i~ p — > 7i Tin amplitudes; the average form factor changes comparatively 
little. The solid curve passes close to the CELSISUS 7r° tt° point [0, but the 
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Figure 2: Total cross sections for the pd — > 3 He(7T7r)° reactions, divided 
by Q 2 as functions of the excess energy Q. The predicted solid and broken 
curves refer to I = 7r° tt° and 7r + 7r~ production respectively, as do the 
closed and open circles from CELSIUS J/J. The triangle is the published 



MOMO 7T+7T- data point 0. The near-threshold IUCF |15j tt+ tt~ point 
(square) is strongly influenced by Coulomb distortion. 



broken one is significantly too low, indicating the presence of some / 



7r + 7r~ production. The IUCF point was obtained at Q = 0.67 MeV [15 1, and 
hence can be assumed to be purely s-wave, though it is heavily influenced by 
Coulomb effects. The comparison of our total cross section predictions with 
the MOMO 7r + 7r~ point at Q = 70 MeV || would suggest that it is mainly 
1 = pairs which are being produced, though this is at variance with the 
strong dependence observed on the angle of the ir + tt~ relative momentum. 
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The 1 = distributions expected at the CELSIUS energy are illus- 
trated in Fig. 3 and these demonstrate the shift to higher masses as compared 
to phase space, which is apparent in the tx N — > nn N input. These experi- 
mental data have insufficient statistics to draw definitive conclusions on 
the shape of the spectrum. It should be noted that the CELSIUS integrated 
cross section points shown in Fig. 2 are mainly determined by the higher 
statistics of their inclusive measurement, which was carried out simultane- 
ously . On the other hand, a low m T7r suppression in the 70 MeV MOMO 
data is clear in their high statistics exclusive 7r + 7r~ production results shown 
in Fig. 4. 

The MOMO data |J show a strong dependence upon the angle 9 n7Tp be- 
tween the relative tttt momentum and that of the beam direction. Taken 
together with the suppression of events at low m^, this suggests the produc- 
tion of / = 1, i = 1 7T7T pairs with spin projection m = ±1 along the beam 
direction. Such an interpretation is backed by the group's preliminary data 
on the inclusive pd — > 3 HeX° reaction ||, which indicate a tt° 7r° production 
rate less than half of that predicted in Fig. 2. 

The production of I = 0, £ = nn pairs would give no dependence upon 
Ottttp, though an interference with an / = 0, £ = 2 contribution could lead 
to such a variation. However, there is no sign of any effect of this kind in 



7i p — > 7r + 7T n |L2|. Since tt + tt p- waves are so small in 7r p — > 7r + 7r n 



near threshold [fL2]| , any simple extension of our model to include p-wave 
production cannot lead to tttt p-wave dominance. 
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Figure 3: Predicted I = tttt effective mass distributions for the pd — > 
3 He 7r + 7r~/7T° 7T° reactions at an incident energy of 477 MeV compared with 
the CELSIUS experimental data J7|. 

For two-pion production near threshold, the intermediate pion in Fig. 1 
gets closer to its mass shell when the 3 He emerges along the direction of 
the initial proton beam and this increases the magnitude of the average 
form factor. The two-step model therefore predicts that, for low m^, the 
dipion should be produced preferentially in the backward hemisphere. This 
effect will, of course, disappear at high masses because the situation then 
approaches one of near-threshold kinematics. 

We have shown that the gross features found in the production of isoscalar 
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Figure 4: Predicted I = 7r + 7r~ effective mass distributions for the pd —>■ 
3 He 7r + 7r~ reactions at an incident energy of 546 MeV, compared with the 
COSY experimental data |J. 

pion pairs in the pd — * 3 He (-7T 7r)° reaction near threshold can be understood 
in terms of the creation of an intermediate virtual pion beam, which in turn 
produces a second meson. The only way that such a model could generate an 
ABC peak in the Q ~ 250 MeV region is if this were already present in the 
tt~ p — > 7t° 7T n input. The parameterisation of the results of the Valencia 
model in Eq. @ corresponds to a parabola in m W7r , whose minimum moves 
to higher values as Q increases. This is due to the enhanced importance 
of the Roper contribution and may leave space at low masses for an ABC 
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effect at higher Q. The question may soon be resolved, because data on 
7i~ p — > 7T° 7T° n at = 750 MeV/c are currently being analysed by the 
Crystal Ball collaboration [16]. 

Further theoretical work is needed to include a more detailed descrip- 
tion of the 7r iV — > 7i 7i N input, though any improvement will, inevitably, be 
rather model-dependent. 
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